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Abstract 


Some  extremely  large  oceanic  trade-wind  cumulonimbus  clouds  extending 
upwards  of  U0,000  ft.  into  a region  of  strong  winds  and  intense  vertical 
shear  have  been  studied  by  means  of  time-lapse  photography.  A simultaneous 
still  picture  of  the  clouds  taken  a known  distance  and  direction  away  from 
the  motion  pictures  permits  triangulation  upon  the  clouds  and  calculation  of 
the  vertical  and  horizontal  motions  of  several  of  the  individual  towers.  By 
means  of  a nearby  radiosonde  observation,  it  is  established  that  the  air 
forming  the  strongest  of  these  towers  must  have  risen  from  near  cloud  base 
with  little  or  no  dilution  of  its  buoyancy  by  mixing  with  the  clear-air  sur- 
roundings . 

The  model  of  Malkus  and  Scorer  (195U)  concerning  the  erosion  of  cumulus 
towers  is  reviewed  and  tested  upon  these  towers  with  satisfactory  results. 

A minimum  horizontal  cloud  dimension  is  apparently  required  for  the  produc- 
tion of  undiluted  towers..  The  horizontal  motion  of  the  clouds  relative  to 
the  air  is  also  estimated  from  the  model  and  tested  by  the  observations  and 
the  net  upward  transport  of  latent-  heat  in  water  vapor  is  calculated  approxi- 
mately. 


I 


INTRODUCTION 


During  the  spring  of  1953,  a joint  observational  program  was  undertaken 
by  the  meteorological  group  from  the  Woods  Hole  Oceanographic  Institution  and 
the  Department  of  Meteorology  of  the  Imperial  College,  London.  The  over-all 
purpose  of  the  program  was  the  investigation  of  the  structure  of  the  low-level 
trade  current  and  the  transports  therein  of  heat,  water,  and  momentum,  with 
considerable  emphasis  upon  continuing  the  study  of  trade-wind  clouds.  Two 
main  observation  points  were  used,  the  first  being  the  snail,  f la h We3t  Indian 
island  of  Anegada  (lB^O'N,  6l*°20,W:  dimensions  two  by  ten  miles)  and  the 
second  a PBY-^A  aircraft,  equipped  as  a meteorological  tool,  which  was  flown 
in  the  near  vicinity.  The  present  paper  reports  some  results  of  the  cloud 
observation  program. 

The  cloud  program  was  divided  into  two  main  parts:  the  first  consisted 

in  airplane  traverses  through  a single  cloud  following  the  methods  developed 
by  the  Woods  Hole  group  in  the  San  Juan  expedition  of  1952  (reported  upon  in 
part  by  Malkus,  195U) ; and  the  second  centered  upon  time-lapse  motion  pictures 
made  from  the  island  of  Anegada.  Two  main  objectives  of  the  photography  were 
to  trace  individual  clouds  through  a large  portion  of  their  life  history  and 
to  provide  information  on  the  interaction  of  groups  of  clouds,  neither  of 
which  could  be  achieved  directly  by  airplane  measurements  alone. 

During  the  observation  period  from  Anegada  a particularly  interesting 
synoptic  situation  was  encountered  on  April  1,  1953.  On  this  occasion  sev- 
eral large  cumulonimbus  clouds  reaching  up  to  well  over  U0,000  ft.  were  ob- 
served far  out  over  the  ocean.  Fortunately,  not  only  were  these  clouds  re- 
corded by  the  time-lapse  camera,  but  sufficient  supplementary  information  was 
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available  to  permit  numerous  calculations  from  the  films.  After  a brief  de- 
scription of  the  equipment  and  techniques  used,  the  maior  portion  of  this 
paper  will  be  concerned  with  these  calculations  and  their  interpretation. 

The  synoptic  situation  was  dominated  by  the  eastward  passage  on  the  pre- 
ceding day  of  a strong  polar  trough.  In  the  Caribbean  area,  the  low-level 
easterlies  were  overlaid  by  westerlies  beginning  at  15,000  ft.  elevation. 
These  westerlies  were  concentrated  in  a shallow,  high-tropospheric  layer 
with  maximum  speeds  of  100  mph  at  an  elevation  of  U0,000  ft.  The  sequence 
of  events  is  illustrated  in  the  time  section  for  San  Jvan,  Puerto  Rico  (120 
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vation,  1220-13U5  LST  (1620- 17U5  OCT)  the  region  was  dominated  by  the  north- 
westerly flow  to  the  rear  of  the  trough.  By  this  time  the  overcast  associ- 
ated with  the  disturbance  had  broken  sufficiently  for  several  huge  cumulo- 
nimbus clouds  to  be  observed  end  photographed  from  Anegada.  Fortunately 
the  same  clouds  were  simultaneously  photographed  with  a still  camera  from 
the  FBI,  which  was  a known  distance  and  direction  from  Anegada.  Trlanv ela- 


tion therefore  permitted  calculation  of  the  distances,  heights,  and  sizes 
of  the  clouds  and  their  alteration  with  time  from  the  motion  picture  film. 
The  resulting  horizontal  and  vertical  motions  could  then  be  compared  with 
the  high-level  wind  and  radiosonde  observations  which  were  taken  a little 
over  an  hour  earlier  by  the  Weather  Bureau  at  San  Juan. 


II . METHODS  OF  OBSERVATION  AND  CALCULATION 

The  time-lapse  motion  picture  camera  on  Anegada  was  a Bolex  16  mm, 
equipped  with  an  extreme  wide-angle  (58.3°  horizontal,  lUl.3°  vertical)  lens. 
Exposures  were  made  automatically  at  the  rate  of  one  frame  every  3.6  seconds 
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or  1000  frames  per  hour.  The  camera  was  mounted  on  a tripod  and  aimed  toward 
22°  from  true  North  during  the  entire  run. 

The  still  camera  ip  the  FBY  was  a conventional  Speed  Graphic  (50.0°  hori- 
zontal, U2.0°  vertical)  which  was  mounted  on  gimbals  so  that  it  always  re- 
mained level  and  was  pointed  at  right  angles  to  the  fuselage.  At  the  time 
the  still  pictures  of  the  clouds  were  taken  (1325  LST)  the  aircraft  was  fly- 
ing at  6000  ft.  elevation  toward  25*3°  True  and  was  making  a speed  of  approxi- 
mately 130  mph  over  the  ground.  Since  it  was  observed  over  Anegada  on  this 
heading  at  13lU  or  11  minutes  earlier,  a base  leg  for  triangulation  purposes 


of  2h  miles  is  obtained.  This  distance  was  checked  by  additional  knowledge 


of  the  time  when  the  aircraft  reached  a location  due  north  of  San  Juan.  The 
base  leg  derived  thus  by  interpolation  came  out  23  miles  or  within  5%.  Cor- 
rection of  the  base  lag  for  wind  drift  of  the  aircraft  was  made,  but  was 
negligibly  small. 

Figure  2 shows  the  Speed  Graphic  (PBY)  photograph  used  in  the  triangula- 
tion. The  time-lapse  picture*  vere  printed  at  intervals  of  1.2  minutes,  and 
the  prints  numbered  from  the  beginning  to  the  end  of  the  run  (1-58).  Se- 
lected prints  from  this  series  are  reproduced  in  Figures  3 and  U.  Print  UU 
was  picked  as  the  frame  coincident  with  Figure  2.  Triangulation  was  first 
performed  fox'  the  three  points  marked  with  inverted  arrows  (on  Figure  2 and 
Print  Uii)  on  the  far  cloud,  hereafter  to  be  called  Cloud  I.  The  angle  of 


each  point,  from  the  optical  axis  of  the  camera  was  found  by  using  the  pro- 
portional distance  of  the  picked  point  from  the  center  to  the  edge  of  the 
negative.  The  half-width  of  the  negative  and  focal  length  of  the  lens  were 


accurately  known  and  these  gave  the  tangent  of  the  total  half-angle  covered 
by  the  camera. 
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Results  of  this  and  a similar  triangulation  on  one  point  on  the  near 
cloud  (called  Cloud  II)  are  shown  in  Figure  5.  By  laying  out  t~o  angles  for 
each  point  at  the  ends  of  the  base  leg  on  a scaled  map,  it  was  determined 
that  the  three  points  in  Cloud  I lay  not  only  in  a straight  line  but  almost 
exactly  parallel  to  the  image  plane  of  the  motion  picture  camera  (i.e.  in 
the  plane  of  the  photographs)  and  a perpendicular  distance  from  it  of  78 
miles.  The  point  on  Cloud  II  was  a perpendicular  distance  of  66  miles  from 


the  image  plane.  Aftsx  discussing  the  pOoSible  errors  in  triangulation,  we 
may  carry  cut  calculations  from  this  point  forward  using  the  motion  picture 
frames  alone* 


In  the  case  of  Cloud  I,  for  which  all  the  significant  calculations  were 
performed,  the  image  location  never  exceeded  an  angle  of  + 10°  from  the  opti- 
cal axis,  hence  errors  in  distances  and  separations  due  to  distortion  qy  the 
lens  may  be  totally  neglected  compared  to  errors  in  triangulation.  An  error 
of  2.2  miles  ir.  the  base  line  (9^)  would  give  an  8^  error  in  the  distance  to 
Cloud  I,  and  it  is  doubtful  that  the  error  in  the  base  line  « as  this  great. 
The  largest  source  of  error  lies  in  the  determination  of  the  angles.  A 2° 
error  in  any  one  of  the  angles  would  produce  a 7 mile  or  9%  error  in  dis- 
tance to  Cloud  I.  All  the  angles  save  one  may  be  shown  reliable  within  this 
accuracy.  The  greatest  uncertainty  lies  in  the  angle  between  the  optical 
aiis  of  the  Speed  Graphic  and  the  heading  of  the  PBT,  which  was  adjusted  to 
90°  by  sighting  along  the  trailing  edge  of  the  wing  ( fortunately  constructed 
at  exactly  90°  to  the  aircraft  axis).  That  the  errors  in  angle  were  very 
small,  however,  is  suggested  by  several  independent  checks  upon  the  triangu- 
lation . The  first  is  that  all  three  points  on  Cloud  I lie  in  a straight 
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line  along  the  plane  290°-110° « as  an  analysis  using  the  wind  profile  shows 
must  be  the  case.  The  second  is  that  triangulation  gives  horizontal  separa- 
tions on  the  map  between  these  points  departing  by  much  less  than  10%  from 
the  separations  measured  on  the  photographs,  using  the  derived  distance  to 
the  camera  of  78  miles.  The  final  eheck  is  that  triangulation  on  Cloud  II 
gives  its  heights  and  velocities  as  extremely  close  to  those  of  comparable 
portions  of  Cloud  I. 

The  vertical  angular  distance  (on  Print  14+)  between  the  horizon  and  the 
farthest  left  point  on  Cloud  I is  determinable  from  the  ratio  of  this  dis- 
tance to  the  vertical  half -width  of  the  negative  and  knowledge  of  the  focal 
length  of  the  lens.  Using  the  horizontal  distance  of  78  miles  and  the  tan- 
gent of  this  angle,  the  vertical  elevation  of  this  point  above  the  horizon 
was  found  to  be  27,000  ft.  or  8.21+  km.  This  length  could  then  be  used  as  a 
basic  distance  scale  for  the  calculations  to  follow.  In  order  to  obtain  the 
elevation  cf  this  point  above  sea  level,  however,  a correction  had  to  be 
xm-ruauceu  uuc  l»u  one  oui'vdtui'6  ui  uxc  earth.  For  a distance  of  78  miles, 
this  proved  to  be  3+000  ft.  or  1.2  km.  In  presenting  actual  altitudes  of 
all  points  on  Cloud  I,  this  correction  has  been  added. 

In  order  to  apply  the  distance  scale  obtained  from  Print  1+1+  to  other 
parts  of  the  cloud  and  to  motions  of  portions  of  the  cloud  towers  throughout 
the  entire  70-minute  run,  it  was  first  necessary  to  ascertain  whether  the 
cloud  had  any  component  of  motion  into  or  out  of  the  plane  of  the  photograph 
and  whether  it  might  develop  any  slope  normal  to  that  plane.  Inspection  of 
the  San  Juan  radio-wind  observation  (see  Figures  1 and  6)  snows  that  the 
wind  directions  throughout  the  height  range  of  the  calculations  (30,000- 
1+7,000  ft.)  departed  rarely,  and  never  more  than  by  10°,  from  the  plane  of 


the  photographs.  Although  the  wind  departed  significantly  from  this  direc- 
tion at  lower  levels,  it  may  be  shown  by  consideration  of  even  the  minimum 
horizontal  drag  forces  ordinarily  acting  on  cumulus  that  the  component  of 
cloud  motion  normal  to  the  photographs  should  have  become  vanishingly  small 
by  9.5  'em.  Thus  the  distance  scale  for  Print  1*1*  was  used  throughout  without 
correction,  with  extremely  consistent  results.  The  Appendix  considers  the 
applicability  of  the  San  Juan  soundings  to  the  time  and  location  of  the  clouds, 
and  shows  that  possible  changes  in  wind  and  thermal  structure  were  probably 
too  small  to  disturb  the  calculations. 

III.  RESULTS  OF  SOME  VELOCITY  CALCULATIONS 

Some  of  the  basic  dimensions  of  Cloud  I are  readily  obtainable  from 
Print  1*1*.  At  31,000  ft.  (9.5  km),  the  width  of  the  left-hand  tower  (contain- 
ing the  left-most  point  triangulated)  is  20,500  ft.  or  6.3  km.  The  great- 
est elevation  above  saa  level  attained  by  the  cloud  is  reached  by  the  streamer 
extending  out  to  the  right  which  has  a maximum  height  of  about  15.3  km  or 
1*6,600  ft.  and  a horizontal  extent  of  1*5  miles.  The  striking  feature  of 
this  cloud’s  changes  with  time  on  Prints  1-58  (70  minutes)  is  its  production 
of  successive  cowers,  each  of  which  is  seen  to  rise  vertically,  spread  lat- 
erally into  a long  streamer  and  then  dissipate  gradually.  On  four  of  the 
more  vigorous  of  these  towers,  vertical  and  horizontal  speeds  were  calculated 
on  the  topmost  point  as  it  rose.  The  calculation  for  the  tower  marked  X on 
Prints  37,  1*0,  and  1*1*  (Figures  3d,  l*a,  ard  l*b,  respectively)  is  probably 
the  most  accurate,  due  first  to  coincidence  with  time  of  triangulation  and 
second  to  the  ease  with  which  a readily  identifiable  portion  of  the  tower 
could  be  traced  from  one  frame  to  the  next.  The  result  of  the  calculation 
is  presented  in  Table  I. 


- 8 - 


Tabls  I 

Heights  and  Velocities  for  Tower  X,  Cloud  I 


Print 

Time  Interval 

Height 

Height 

w 

uc 

UE 

u = UE-Uc 

number 

minutes 

ft 

km 

Dips 

mps 

mps 

mps 

33 

2.U 

31,000 

9.U5 

10.5 

17.0 

32 

15.0 

35 

2.U 

35,900 

10.95 

7.1 

28.5 

38.5 

10.0 

37 

■so  onn 

y 4.  vv 

11.96 

2.U 

2.7 

19.5 

39 

3.6 

140,500 

12.35 

1.1 

27.5 

1 n 

U.8 

Ul,200 

12,58 

1.1 

29.0 

U6 

U2,l*00 

12.90 

The 

vertical  rate  of 

ascent,  w, 

of  the 

tower  top 

was  obtained  merely  by 

subtracting  the  two  heights  and  dividing  by  the  time  interval,  and  is  thus  an 
average  between  the  frames  indicated.  The  horizontal  velocity,  uc,  of  the 
tower  top  was  obtained  similarly  by  subtraction  of  two  ouccessivc  horizo  ntal 
positions.  The  external  wind,  ug,  was  read  from  Figure  6 at  a level  midway 
between  the  two  elevations,  at  which  heights  X's  were  plotted  showing  the  cal- 
culated value  of  lie.  The  drop  in  Uq  between  11.96  and  12.35  km  may  be  in- 
ferred to  be  due  to  evaporation  of  the  cloud  streamer  from  its  forward  end. 
From  the  film  and  the  computed  velocities  for  9.50-11.96  km,  however,  it  ap- 
pears that  the  tower  is  moving  without  evaporation  and  hence  its  speed  may- 
be inferred  to  be  the  same  as  that  of  the  updraft  for  the  lowest  two  X*s  on 
Figure  6.  The  relative  horizontal  velocities,  u = ug  - Uq,  are  to  be  used 
later  in  examining  the  frictional  interaction  between  cloud  tower  and  sur- 
roundings. 


Two  additional  cloud  towers  on  Cloud  I proved  traceable  throughout  their 
life  histories.  The  first  of  these  appeared  in  Prints  6-27  and  is  marked  by 
horizontal  bars  on  Figures  3a,  b,  and  c.  This  tower  reached  the  greatest 
heignt,  as  will  be  brought  out  in  the  following  section,  and  exhibited  the 
largest  cross  section  of  all  of  them.  The  pertinent  calculations  are  pre- 
sented in  Table  II.  Due  to  the  considerable  width  of  the  tower,  its  hori- 
zontal motion  could  not  be  determined  with  great  reliability  after  Print  9. 

Table  II 

Heights  and  Velocities  for  Tower  I,  Cloud  I 


Print 

Time  Interval 

Height 

Height 

w 

uc 

number 

minutes 

ft 

km 

mps 

mps 

6 

3.6 

*■*  > OAr\ 

J4, ouu 

T A Cf\ 

_LU. 

8.1* 

2$. 2 

9 

3.6 

1 r\  *7A/-N 

mi.  1 UU 

12.1*0 

1.9 

28.2 

12 

3.6 

1*2,000 

±2.80 

1.1* 

3$ 

1$ 

3.6 

1*3,100 

13.10 

1.2 

32 

18 

10.8 

1*3,900 

13.3$ 

0.3 

indet  erminate 

27 

1*1*,  1*00 

13.$$ 

These  values  of  uc  are  shown  as  horizontal  bars  on  Figure  6 at  the  mean 
of  the  two  indicated  elevations. 

The  third  tower  was  considerably  weaker  than  the  two  previous  ones,  and 
dissipated  at  a considerably  lower  elevation.  It  appeared  on  Prints  21-3$ 
and  is  shown  by  the  inverted  triangle  on  Figure  3c.  The  pertinent  calcula- 
tions are  presented  in  Table  III. 
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Table  III 

Heights  and  Velocities  for  Tower  Q,  Cloud  I 


Print 

'fime  Interval 

Height 

Height 

w 

uc 

UE  u 

= U£-uc 

number 

minutes 

ft 

km 

mps 

mps 

mps 

mps 

21 

2.4 

31,200 

9.5 

7.0 

13.5 

22.5 

9.0 

23 

2.4 

34,500 

10.5 

3.5 

24.0 

37.0 

13.0 

25 

2.4 

36,100 

11.0 

2.1 

24.0 

39.0 

15.0 

27 

2.4 

37,200 

11.32 

1.1 

26.0 

OO 

«-✓ 

7.2 

37,650 

11.49 

0.6 

20.0 

35 

38,850 

11.84 

The 

values  of  are  plotted 

at  the  appropriate  heights 

on  Figure 

6 using 

inverted  triangles. 

In  addition  an  early  tower,  called  W,  on  Cloud  I was  studied,  but  it  was 
closing  its  life  cycle  at  the  beginning  of  the  run  and  hence  it  is  not  pre- 
sented here  in  detail.  Its  horizontal  velocities,  however,  are  entered  as 
squares  in  Figure  6,  and  it  is  shown  by  the  streamer  with  a square  above  on 
Figures  3a  and  b.  The  streamer  on  Cloud  II  was  also  followed  (shown  by  tri- 
angles on  Figures  3c,  d,  and  4a,  b,  c,  and  d)  and  its  horizontal  velocities 
are  given  by  the  triangles  at  the  appropriate  heights  on  Figure  6.  This 
calculation  was  continued  up  to  higher  elevations  than  were  the  towers  on 
Cloud  I merely  because  the  cloud  fragment  was  easier  to  detect.  Some  of  the 
streamers  from  Cloud  I attained  equally  great  elevations. 

Now  a physical  interpretation  of  some  of  the  foregoing  observations  may 
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be  possible.  It  is  important  to  note  that  the  three  towers  described  here 
underwent  an  initial  period  or  rapid  rise  which  petered  out  into  a period  of 
much  slower  rise  somewhere  between  about  11.0  and  12.5  km.  In  the  case  of 
the  most  vigorous  towers,  X and  Y,  close  inspection  of  the  films  reveals 
that  this  transition  took  place  very  near  to  12.2  and  12.6  km,  respectively, 
where  the  appearance  of  the  towers  also  changed  markedly  and  sharp  outlines 
became  lost.  This  transition  could  be  in  no  way  related  to  the  transition 
between  super-cooled  water  and  ice  crystals  which  must  have  occurred  at  no 
higher  than  10=9  km,  where  even  the  most  buoyant  possible  cloud  element 
reached  -Ii0.0°0 . 

The  San  Juan  radiosonde  observation  for  1500  GOT  is  reproduced  in  Fig- 
ure 7-  The  sounding  has  a most  striking  feature.  It  shows  that  any  cloud 
element  which  is  finally  to  reach  the  vicinity  of  12.5  km  with  vanishing 
(and  not  negative)  buoyancy  must  have  risen  moist  adiabatically  and  undi- 
luted all  the  way  from  cloud  base  at  950  mb.  Any  mixing  or  dilution  by  the 
clear-air  surroundings  on  the  way  up,  sv  . a small  fraction  of  that  encoun- 
tered by  Stommel  (19U7),  Malkus  (1952,  195k)  and  others  in  smaller  trade 
cumulus  would  result  in  prohibitively  large  negative  buoyancies  at  this 
level.  The  next  section  is  intended  to  shed  light  upon  this  point. 

IV.  THE  EROSION  AND  LIFE  CYCLE  OF  THE  CLOUD  TOWERS 

a)  The  Model  To  Be  Used 

The  successive  production  by  a cumulus  of  isolated  towers,  and  their 
subsequent  dissipation,  is  a situation  clearly  inappropriate  for  the  appli- 
cation of  a steady-state?  cloud  model.  A theoretical  framework  for  describ- 
ing the  erosion  of  such  transient  towers  has,  however,  recently  been  worked 
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out  by  Malkus  and  Scorer  (195U)  in  a paper  hereinafter  referred  to  as  I.  The 
present  observations  provide  material  for  the  further  testing  and  possible 
extension  of  their  model.  These  authors  treated  cumulus  towers  as  isolated, 
buoyant  air  bubbles  which  eroded  as  they  rose,  each  closing  its  individual 
life  cycle  when  totally  eaten  away. 

The  study  of  large  buoyant  bubbles  in  application  to  geophysical  problems 
was  commenced  by  G.  I.  Taylor  (Davies  ana  Taylor,  1950 j Taylor,  1950).  Davies 
and  T ay lor  studied  air  bubbles  rising  through  water  and  showed  that  the  ob- 
served limiting  velocity  of  rise  could  be  obtained  by  assuming  potential  flow 
around  the  spherical  upper  surface  of  the  bubble.  This  limiting  velocity  is 
reached  because  the  buoyancy  of  the  bubble  is  balanced  by  the  drag  exerted  by 
the  surrounding  fluid  as  it  is  deflected  around  the  cap.  Taylor  (1950)  ap- 
plied his  results  to  buoyant  air  bubbles  in  air  in  the  case  of  the  atomic 
bomb  clcud  with  surprisingly  high  accuracy.  Scorer  and  Ludlam  (1953)  sug- 
gested that  the  bubble  model  could  be  applied  with  value  to  ordinary  convec- 
tive elements,  such  as  tiiSiiuulS  3 nd  cumulus  clouds,  and  presented  evidence 
that  cumulus  towers  showed  both  the  required  rounded  tops  and  approximated 
a limiting  velocity.  They  added  the  idea  of  erosion,  hypothesizing  that  air 
bubbles  in  air  would  retain  the  near-spherical  cap  by  a continuous  shedding 
of  their  outer  skin,  which  sank  relatively  and  became  incorporated  into  the 
wake.  That  spherical- capped,  bubble-like  buoyant  elements  may  arise  sponta- 
neously in  a convecting  fluid  nas  been  noted  by  W.  Malkus^  in  his  extensive 
laboratory  experiments  (19b4)  on  a fluid  heated  from  below* 

In  paper  I,  a quantitative  law  for  the  arosion  of  atmospheric  bubbles 
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Conversations  with  the  writers. 
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was  put  forward  and  tested.  This  law  states  that 

dR3  2 

= -3ER  gB  (1) 

dt 

where  R is  the  radius  of  curvature  of  the  bubble,  E is  a constant  of  propor- 
tionality called  the  erosion  constant,  and  gB  is  the  buoyancy  acceleration, 
where  B = ZX  /°  I i0  , ZX  Z-3  being  the  bubble ' s density  deficiency  with  respect 
to  that  of  the  environment  of  density  z0  . 

If  the  buoyancy  of  the  bubble  is  constant,  that  is,  when  the  environ- 
ment has  a wet  adiabatic  lapse  rate,  we  may  readily  integrate  (1).  Using 
tne  boundary  condition  R = 0 at  t = 0,  and  considering  negative  time  decreas- 
ing during  the  bubble's  life,  this  integration  gives 


R = -EgBt  (2) 

or  a linear  decrease  of  radius  with  time  is  hypothesized.  This  erosion  law 
was  combined  with  an  equation  of  motion  for  unit  mass  of  the  bubble,  namely 


dw 


+ Kw  = gB 


(3) 


where 


r» 

K = — (3a) 

UH 


A differential  equation  of  this  form  may  be  shown  applicable  to  cloud  towers 


provided  only  that  they  retain  a quasi-spherical  top  and  do  not  erode  so  fast 


that  the  actual  velocity  departs  radically  from  the  limiting  velocity  wQ. 
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The  limiting  Telocity  may  be  found  from  equation  (3)  by  setting  dw/dt,  here- 
after called  equal  to  aero,  via. 

2 1 

wQ  = - (gBR)2  (U) 


whien  is  the  same  expression  derived  and  tested  by  Davies  and  Taylor  (193>0) 
when  B 1. 

When  the  erosion  law  (2)  is  substituted  into  the  differential  equation 
(3)  and  the  result  solved  for  E we  have 


3 X (-t)i  gB  (1  - w/gB)^ 

- E 2 = G = 

2 w 


(5) 


If  the  hypo the siaed  laws  are 


valid,  3/2 


or  G should  be  constant  for 


XI 

ou-l  j.  j-oxoiv-va-j 


iiivxL  v xuucLL 


bubbles.  Equation  (3>)  is  ideally  suited  to  test  by  calculations  from  time- 
lapse  motion  pictures,  if  the  buoyancy  gB  is  known  or  can  be  estimated  and 
does  not  vary  too  rapidly  with  time.  It  was  shown  in  I (see  Section  lib  and 
Appendix  I,  loc.  cit.)  that  equation  (5)  could  be  used  with  varying  values  of 
buoyancy  provided  that  the  buoyancy  did  not  vary  by  more  than  1/3  ir  100 
second?  * An  inc!*6H5?  buoysncy  at  this  ovsr^stinists  0?  s.s 

much  as  20%.  It  should  be  pointed  out  that  from  the  level  of  their  forma- 
tion upward,  the  cloud  bubbles  are  hypothesized  to  be  rising  along  a wet 
adiabatic,  mixing  with  the  environment  being  achieved  by  erosion  and  wake 


production,  the  core  remaining  undiluted  until  its  final  disappearance.  The 
bubbles  studied  in  I were  found  to  originate  at  different  levels,  some  over 
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land  forming  at  or  near  the  ground,  others  at  cloud  base,  and  still  others 
up  to  the  highest  level  within  the  cloud  where  the  external  lapse-rate  ex- 
ceeded moist  adiabatic.  No  trade  cumulus  bubbles  were  found  in  that  study 
which  originated  below  cloud  base.  The  ones  formed  at  higher  levels  were 
generally  composed  of  smaller  ones  plus  wake  material,  and  thus  were  more 
"dilute"  or  had  less  buoyancy  than  those  forming  lower  down. 

Tests  performed  in  I upon  small  and  medium-sized  cumulus  bubbles  showed 
that  G was  constant  to  an  excellent  degree  of  accuracy,  both  during  the  life 
cf  a single  cloud  tower  and  from  tower  to  tower,  regardless  of  location  or 
external  conditions.  The  mean  value  of  G came  out  0.2,  giving  an  E of  about 
5>0  seconds.  The  highest-penetrating  bubble  studied  was  a ix-ade  cumulus 
tower  reaching  16,000  ft.  which  did  not,  however,  become  glaciated. 


b)  The  Life  Cycles  of  Individual  Towers 

The  cloud  towers  in  the  present  study  not  only  reached  great  heights, 
but  penetrated  tremendously  high  wind  shears  and  definitely  contained  ice 
crystals  throughout  a significant  portion  cf  their  life  histories _ It 
should  be  of  considerable  importance  to  determine  whether  equation  (5)  and 
the  constant  value  of  E obtained  for  the  bubbles  studied  in  paper  I give 
here  also  a consistent  relation  between  velocity,  acceleration,  buoyancy, 
and  time.  If  bu,  it  will  be  of  interest  to  see  how  the  present  observations 
may  be  interpreted  and  clarified  by  means  of  this  model.  Most  of  the  pres- 
ent towers  have  the  advantage  of  being  well  isolated  from  interference  by 
one  another  due  to  great  wind  shear  and  thus  form  an  especially  valuable 
oest  case. 


To  make  this  test,  the  original  motion  picture  film  was  projected  on  a 


- 16  - 


microfilm  viewer,  from  which  the  heights  of  the  bubble  tops  were  read  off  at 
every  frame  (3.6  sec).  Although  the  vertical  resolution  was  still  smaller 
than  desirable,  there  being  commonly  only  about  1 cm  height  difference  be* 
tween  the  emergence  of  a bubble  and  its  death  as  an  entity,  the  results  of 
drawing  a smooth  curve  through  the  points  were  extremely  satisfactory.  The 
measurements  and  method  of  calculation  of  velocities  and  accelerations  for 
Tower  X are  shown  in  Figure  8.  A dual  criterion  was  used  for  picking  the 
time  of  complete  erosion,  when  (-t)  was  to  be  set  equal  to  zero,  as  described 
earlier  in  I.  This  time  was  chosen  when  the  levelling  off  of  the  ascent  was 
indicated  on  the  graph,  coincident  with  a fuzzy  appearance  and  loss  of  clear- 
cut  edges  showing  itself  on  the  film.  In  applying  and  testing  equation  (£) 
on  Tower  X,  therefore,  the  quantities  w,  w and  -t  were  taken  from  the  curves 
shown  in  Figure  8.  It  was  then  desired  to  see  whether  reasonable  and  ob- 
servationally  supportable  values  of  gB,  the  bubble 1 s buoyancy,  could  give 
rise  to  values  of  G which  both  remained  nearly  constant  throughout  the  bub- 
ble’s life  and  averaged  about  0.2  sec-^  as  established  for  the  bubbles  in 
paper  I.  When  the  values  of  buoyancy  are  chosen  with  the  aid  of  the  sound- 
ing in  Figure  7,  under  assumptions  to  be  described  presently,  the  resulting 

calculations  of  G are  obtained  which  are  shown  in  Table  IV.  Under  these 

_ 

conditions  the  average  G,  or  G,  is  0,19$  sec  2 and  the  percentage  spread 
in  G,  called  , is  only  2 $%, 

The  values  of  buoyancy  in  the  column  marked  gB  in  Table  IV  may  be  ob- 
tained for  this  bubble  from  the  sounding  in  Figure  7 under  several  assump- 
tions, the  most  crucial  being  that  the  San  Juan  radiosonde  observation  is 
applicable  to  the  clear-air  environment  of  the  clouds . This  point  is 
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Table  17 

Tower  I,  April  1,  1953 
Height  range  studied  10.U2-12.2U  km 


Observa- 

tion 

-t 

(-t)* 

w 

• 

w 

gB 

w/gB 

(1-w/gB)* 

r\ 

U 

number 

sec 

sec2 

mps 

cm/sec^ 

cm/sec^ 

sec  2 

1 

288 

17.0 

11.7 

-2.2 

1U.0 

-.157 

1.08 

0.22 

2 

270 

16.5 

11.2 

-2.2 

13.3 

-.165 

1.08 

0.21 

3 

252 

15.9 

10.5 

-2.2 

12.7 

-.173 

1.08 

0.21 

U* 

23U 

15-3 

9.8 

-2.2 

11.5 

-.190 

1.09 

0.20 

5 

216 

1U-7 

9.2 

-2.2 

10.U 

-.210 

1.10 

0.18 

6 

198 

1U.1 

8.3 

-2.8 

9.0 

-.310 

1.15 

0.18 

7 

180 

13.5 

7.2 

-2.8 

8.8 

-.320 

1.15 

0.19 

8 

162 

12.8 

6.U 

-3.0 

8.6 

-.350 

1.16 

0.20 

9 

1UU 

12.0 

5.8 

-3.0 

8.U 

-.360 

1.17 

0.20 

10 

126 

11.2 

5.5 

-2.2 

8.2 

-.270 

1.13 

0.19 

11 

108 

10. U 

5.3 

-2.2 

8.0 

-.280 

1.13 

0.18 

12 

90 

9.5 

U.7 

-3.0 

7.8 

-.385 

1.18 

0.19 

13 

72 

8.U 

U.U 

-3.0 

7.6 

-.395 

1.18 

0.17 

lU 

5U 

7.U 

3.9 

-3.8 

7.6 

-.500 

1.23 

0.18 

15 

36 

6.0 

3.0 

t Q 
— U#'-' 

n Cs 
l 

Aon 

— • 

1.28 

0.19 

16 

18 

U.2 

2.0 

-6.2 

7.6 

-.820 

1.35 

0.22 

^Bubble 
as  1^3. 

diameter 
3 "<m 

measured 

G = 

A = 

0.195 

25* 
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established  in  the  Appendix.  The  second  assumption  is  that  the  ascent  curve 
of  the  bubble  forming  Tower  X is  very  closely  given  by  the  dashed  (right-hand) 
moist  adiabat  on  Figure  7,  which  originates  from  the  environment  sounding  at 
cloud  base  level.  That  this  curve  delineates  the  maximum  possible  buoyancy 
for  this  tower  is  supported  by  other  work  on  oceanic  trade  cumuli  (Malkus, 

19  5U)  showing  that  cloud  buoyancies  at  cloud  base  level  rarely  exceed  that- 
contributed  by  water  vapor  alone.  It  will  be  noted  that  this  adiabat  re- 
crosses the  sounding  curve  just  above  12.5  km,  leaving  the  bubble  some  buoy- 
ancy at  its  demise  at  12.21*  km.  Thus  this  adiabat.  according  to  the  bubble 
model,  also  represents  minimum  buoyancy  for  this  tower  since  isolated  bub- 
bles cannot  continue  to  rise  without  positive  buoyancy.  In  proceeding  from 
the  virtual  temperature  differences  between  the  dashed  adiabat  and  the 
sounding  curve  to  the  values  of  gB  presented  in  Table  IV  a negative  correc- 
tion to  the  buoyancy  due  to  the  presence  of  suspended  hydrometeors  has  been 
introduced.  This  correction  allows  for  about  1*.5  gm/m"  of  water  or  ice 
particles  at  10.1*2  km,  diminishing  to  2 gm/m^  at  the  bubble’s  demise  at 
12.21*  km.  This  correction  amounts  to  nearly  30$  of  the  total  buoyancy  and 
is  both  uncertain  and  arbitrary:  it  is  simply  delimited  by  that  reduction 

in  buoyancy  necessary  to  proceed  from  the  dasned  adiabat  to  the  values  of 
gB  in  the  Table  which  produce  values  of  G consistent  with  tuuse  in  paper  I. 
That  these  figures  are  not  unreasonable,  however,  may  be  seen  as  follows: 
if  all  hydroineteors  were  being  retained  in  the  bubbles  from  cloud  base  up- 
ward,  the  water  or  ice  content  at  10.5  km  would  be  about  18  gm/nv  . These 
clouds  were,  however,  observed  to  be  precipitating  heavily  from  their  bases. 
Tims  a hydrometeor  content  from  10-25$  this  figure  is  plausible,  especially 
when  it  falls  within  the  range  of  observed  water  contents  in  large  cumuli. 
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The  latent  heat  of  freezing  has  been  ignored  in  the  present  computation,  the 
assumption  being  that  most  of  the  freezing  took  place  below  10.1*2  km  where 
the  bubbles  began  to  be  tracked. 

Tower  Y was  treated  in  a similar  manner  as  Tower  X and  the  results  are 
presented  in  Table  V.  The  buoyancy  used  in  this  table  was  obtained  also 
using  the  dashed  adiabat  in  Figure  7,  with  the  same  magnitude  correction  for 
weight  of  hydrometeors.  It  should  be  pointed  out  that  Tower  Y was  composed 
of  several  bubbles  of  the  size  of  the  one  tabulated,  which  accounts  for  its 
great  width.  In  Figure  3a,  four  of  these  are  plainly  visible.  The  bubble 
studied  is  the  right-hand  one  which  ascended  first.  However,  it  will  be 
noted  that  the  arrival  of  successors  disturbed  the  constancy  of  G,  especially 
from  t = -210  sec.  on.  This  was  due  to  overly  large  ascent  speeds  acquired 
by  this  bubble  due  to  being  pushed  up  by  successors,  resulting  in  a too  lo« 
value  of  G.  The  same  effect  was  noted  in  several  cases  in  paper  I. 

Next,  one  of  the  weaker  towers,  called  Q,  was  studied  in  the  same  man- 
ner and  the  results  are  presented  in  Table  VI.  The  significant  feature  of 
Table  VI  is  that  a considerably  lower  buoyancy  had  to  be  chosen  for  this 
bubble  to  give  the  same  value  of  G as  the  others.  The  buoyancy  for  this 
bubble  was  taken  from  the  moist  adiabat  shown  by  the  dotted  line  (left-hand 
adiabat)  on  Figure  7,  with  a correction  for  weight  of  hydrometeors  of 
3 gm/m^  at  9.72  km  diminishing  to  2 gm/nP  at  its  maximum  height  of  10. 9U  km. 
The  resulting  values  of  gB  given  in  Table  VI  average  0.99  times  the  values 
used  for  X and  Y in  the  height  range  10.U2-10.9i*  km.  This  latter  bubble 
(q)  possessed  lower  ascent  rates  and  perished  sooner  than  the  preceding  two 
described.  It  can  now  be  seen  that  it  must  have  been  a considerably  more 
dilute  bubble,  its  level  of  origin  being  perhaps  about  9.9  km  or  19,300  ft. 
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Table  V 

Tower  T,  April  1,  1953 
Height  range  studied  10.25-12.60  km 


Observa- 

(-t) 

(-t)* 

w 

w 

gB 

w/gB 

(1-w/gB)2 

G 

tion 

number 

see 

sec^ 

mps 

cm/sec^ 

cm/sec^ 

4 

sec  * 

1 

318 

17.8 

11.7 

0 

15.5 

0 

1.00 

0.23 

2* 

300 

17.  1* 

ll.  1* 

-0.6 

H*.o 

- .ol* 

1.02 

0.22 

3 

282 

16.8 

10.9 

-0.9 

13.3 

- .07 

l.ol* 

0.21 

h 

261* 

16.3 

10.5 

-1.5 

12.7 

- .12 

1.06 

0.21 

5 

2U6 

15.7 

10.3 

-2.2 

11.5 

- .19 

1.09 

0.19 

Z 

ooo 

if'  n 

-» r\  r\ 

-2.8 

1 A A 

O Z 

A TO 

u 

C.C.SJ 

-LU.  V 

- .cO 

J.  * ic 

U#x,7 

7 

210 

ili. 5 

9.8 

-3.1 

9.0 

- .35 

1.16 

0.16 

8 

192 

13.8 

9.2 

-3.1 

8.6 

- .36 

1.17 

0.15 

9 

17U 

13.2 

8.6 

-3.7 

8.1* 

- .1*1* 

1.20 

0.16 

10 

156 

12.5 

7.5 

-3.7 

8.0 

- .1*6 

1.21 

0.16 

11 

138 

11.8 

6.7 

-3.7 

7.6 

- .U9 

1.22 

0.16 

12 

120 

11  = 0 

6.7 

-3.7 

7.6 

- .1*9 

1.22 

0.15 

13 

102 

10.1 

6.1* 

-3.7 

7.6 

- .1*9 

1.22 

0.15 

lii 

81* 

O o 
y • *• 

6.0 

-i*.3 

7.6 

- .56 

1.25 

0.15 

15 

66 

6,1 

5.o 

-5.3 

7.0 

- .76 

1.33 

0.15 

16 

1*8 

6.9 

3.6 

-5.6 

4 c: 

- .90 

1 ofi 

0.17 

17 

30 

5.5 

2.2 

-6.0 

r' 

-1.20 

1.1*9 

0.19 

l8 

12 

o d 

1.0 

-7.0 

l*.o 

-1.75 

1.66 

0.23 

^Bubble  diameter  measured 
as  2.1  km 


G = 0.180 
A = h$% 
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Table  VI 

Tower  Q,  April  1,  1953 


Height 

range 

studied  9. 

72-10. 9U  km 

Observa- 

tion 

-t 

(-t)* 

w 

• 

w 

SB 

w/gB 

(1-w/gB)^ 

G 

number 

sec 

x 

sec2 

mps 

cm/sec*' 

, 2 
cm/ sec 

sec  * 

1* 

216 

lii.7 

8.9 

0 

13.0 

0 

1.00 

0.21 

2 

198 

lU.O 

8.6 

1 

o 

• 

v: 

12.5 

-0.07 

1.0U 

0.21 

3 

ioC 

13. U 

8.U 

-1.2 

12.0 

-0.10 

i.o5 

0.20 

U 

152 

12.3 

8a 

-1.6 

11.1 

-0.1U 

1.07 

0.18 

5 

1UU 

12.0 

7.2 

-1.9 

10.5 

-0.16 

1.09 

0.19 

6 

126 

11.2 

6.U 

-2.2 

9.8 

-0.23 

1.11 

0.19 

7 

108 

10.  U 

5.9 

-2.8 

9.2 

-0.30 

l.lU 

0.19 

8 

90 

9.5 

5.5 

-3.1 

8.8 

-0.35 

ia6 

0.18 

9 

72 

8.5 

5.1 

-3.7 

8.U 

-o.UU 

1.20 

0.17 

10 

5U 

7.U 

U.l 

-U.3 

8.0 

-o.5U 

1.2U 

0.18 

11 

36 

6.0 

3.0 

-5.5 
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rather  than  cloud  base.  It  presumably  formed  from  smaller  bubbles  and  dilute 
wake  material  near  this  level,  while  the  other  two  formed  near  cloud  base. 

Its  concomitantly  smaller  required  temperature  excess  (about  1.3°C  compared 
to  about  2.6°C  for  I and  l)  over  the  surroundings  at  $. 9 km  is  consistent 
with  this  point  of  view. 

It  will  be  seen  by  inspection  of  Figures  3 and  ii  that  Tower  Q rose  on 
the  extreme  left  side  of  Cloud  I (upwind  and  upshear).  In  fact  at  the  time 
it  formed,  this  part  of  the  cloud  was  almost  entirely  separate  from  the 
larger  cloud  mass  to  the  right.  The  average  width  of  its  cloud  mass  was 
only  about  1*  km,  while  Tower  I rose  from  a clouu  utdoD  clu  xcoi^u  y hi.  in  width. 
By  the  time  Tower  1 rose  from  the  right-hand  portion  of  the  cloud,  this  sec- 
tion had  grown  until  it,  too,  had  attained  an  average  width  of  about  9 km. 
These  data  might  thus  suggest  that,  under  given  conditions,  a minimum  hori- 
zontal cloud  dimension  is  required  to  permit  emergent  buoyant  elements  to 
rise  all  the  way  from  cloud  base  without  dilution.  Since  th6  undiluted  ele- 
ments will,  in  general,  compose  the  strongest  and  most  penetrating  towers, 
it  may  well  be  that  for  cumulus  to  attain  a sufficient  vertical  thickness 
to  give  rise  to  precipitation,  a comparable  horizontal  size  must  first-  be 
attained. 

To  investigate  these  points  further,  the  other  weak  towers  put  up  by 
Cloud  I were  studied.  There  were  three  such  towers  detestable  on  the  film 
One  (P)  was  too  feeble  and  short-lived  to  be  treated  quantitatively.  One 
of  the  two  remaining,  called  C,  appeared  on  Print  12  (see  Table  II)  just  as 
the  bubble  studied  in  Tower  I was  closing  its  life  cycle.  Although  it  ap- 
peared from  the  same  large  cloud  mass  as  did  Y,  its  buoyancy  was  the  same 
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as  that  of  the  diluted  tower,  Q.  It  had  an  active  lifetime  of  185  secs,  a 
maximum  updraft  speed  of  7.2  mps  and  a height  range  of  10.L-11.2  km. 

The  final  tower  was  the  very  last  one  appearing  on  the  film  sequence, 
called  R.  In  order  to  have  the  same  value  of  G as  the  other  towers,  it  had 
to  be  slightly  dilute  and  to  follow  a moist  adiabat  about  half-way  between 
the  dashed  and  the  dotted  curve  on  Figure  7.  It  covered  a height  range  from 
10.6-11.85  km,  exhibited  a maximum  updraft  speed  of  7.8  mps,  and  a life 
period  of  275  secs,  just  short  of  those  of  the  strongest  towers.  At  the 
time  of  its  emergence,  the  entire  cloud  mass  of  Cloud  I was  apparently  be- 
ginning to  decay,  and  although  the  bubble  came  from  the  main  portion  of  the 
cloud-  it3  horizontal  dimension  was  then  only  a 1 looic  over  7 kiu<  Tiiizs  it 
is  seen  that  although  the  strength  of  the  emerging  towers  and  the  dimension 
of  the  producing  cloud  mass  in  general  appear  to  have  some  relation,  another 
factor  or  factors  must  also  be  operative,  since  the  largest  cloud  mass  pro- 
duced weak  as  well  as  strong  towers.  A suggestion  concerning  one  of  these 
factors  is  made  in  later  paragraphs. 

The  major  features  of  the  towers  studied  are  summarized  in  Table  711, 
which  also  compares  observed  bubble  diameters  with  the  values  of  their  radii 
of  curvature  calculated  from  the  Malkus-Scorer  model.  In  the  table,  the 
measured  diameter  of  each  bubble  is  compared  to  the  calculated  radius  of 
curvature,  using  both  equation  (2)  and  equation  (u),  for  the  3ame  frame  as 
the  measurement.  The  earlier  work  demonstrated  that  D,  so  measured,  and  R 
came  out  very  nearly  equal.  The  excellent  agreement  in  the  present  cases 
gives  further  support  to  the  applicability  of  the  bubble  model.  The  col- 
umn marked  “max  gives  the  approximate  dimension  of  each  bubble  as  it  first 
emerged  from  the  cloud  top,  calculated  by  extrapolation  from  equation  (2). 
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It  will  be  noted  that  the  strong,  undilute  towers,  X and  T,  had  maximum  radii 
in  excess  of  2 km,  while  the  weaker  ones  just  exceeded  one  kilometer. 

Two  of  the  four  weaker  towers  emerged  from  smaller  cloud  masses  while 
the  other  two  emerged  from  the  largest  (<-^9  km)  ones.  In  the  case  of  the 
latter,  namely  C and  P,  one  noteworthy  feature  is  held  in  common.  In  each 
case,  the  weak  tower  followed  its  predecessor's  emergence  after  an  unusually 
short  interval,  7.2  minutes  in  the  case  of  C and  6 minutes  in  the  case  of  P. 
The  strongest  towers,  J and  X,  were  separated  by  more  than  a half  hour,  while 
X and  R were  separated  by  nearly  20  minutes.  If,  as  suggested  by  writers  on 
the  bubble  model,  the  larger  cloudy  bubbles  are  formed  by  aggregation  of 
auirto-lBr  once  at  low  levels  within  the  cloud,  this  fact  becomes  somewhat 
clarified.  A certain  time  interval  is  apparently  necessary  for  a sufficient 
number  of  small  bubbles  to  accumulate  in  the  low  cloud  levels  in  order  to 
form  a large  vigorous  one  which  can  rise  to  emerge  undiluted  from  cloud  top. 

The  last  four  columns  in  Table  VII  may  now  be  examined.  By  "protected 
range"  is  meant  the  vertical  distance  between  the  level  of  bubble  origin 
and  the  cloud  top.  The  level  of  origin  for  each  bubble  was  estimated 
roughly,  from  Figure  7 in  conjunction  with  the  film.  The  "unprotected 
range"  is  the  height  range  studied,  or  the  height  range  achieved  during  the 
active,  emergent,  isola  ted  life  of  the  bubble.  In  the  case  of  those  bub- 
bles which  aid  not.  follow  their  predecessors  too  closely,  the  protected 
range  and  cloud  mass  width  are  very  comparable.  In  the  case  of  all  bubbles 
studied,  the  unprotected  height  range  and  radius  at  emergence,  Rjuajc,  are 
comparable,  as  noted  for  smaller  bubbles  in  I.  It  appears  that  the  pres- 
ence of  ice  crystals  has  no  noticeable  effects  upon  the  bubbles  during  their 
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active  lives,  although  it  may  cause  slower  dissipation  of  their  wakes  and 
remnants. 

The  last  column  on  the  right  shows  that  the  very  vigorous  undilute  bub- 
bles had  active  lives  of  about  300  secs  or  five  minutes  (X  could  have  been 
followed  this  long  on  the  film)  and  successively  weaker  bubbles  had  succes- 
sively shorter  lives.  Bubble  P probably  existed  independently  for  100  secs 
or  less.  In  each  case,  as  Tables  I-III  show,  the  remnants  of  the  towers 
went  on  rising  at  roughly  1 mps  for  about  10  minutes  more,  achieving  heights 
greater  than  13  km.  Then  dissipation  set  in  and  occupied  another  10-1$  min- 
utes without  much  further  ascent  taking  place. 

Since  the  present  discussion  has  in  part  concerned  undiluted  cloud  ele- 
ments, some  rising  all  the  way  from  cloud  base  moist  adiabatically,  it  must 
be  emphasized  that  any  resemblance  to  the  hypotheses  and  results  of  the  par- 
cel method  is  strictly  superficial.  The  parcel  method  concerns  cloud  ele- 
ments which  in  no  way  interact  with  the  surrounding  air  and  which  would  con- 
tinue to  accelerate  upward  wherever  the  parcel  density  is  less  than  that  of 
the  envi-ronmeut.  In  the  present  case,  a rough  calculation  shows  that  par- 
cels becoming  saturated  at  cloud  base  would  emerge  from  cloud  top  at  10  km 
with  upward  velocities  well  in  excess  of  60  mps  and  still  increasing!  Since 
the  bubble,  on  the  other  hand,  is  hypothesized  to  be  travelling  at  very 
neaily  its  limiting  velocity  at  all  times  fat  least  after  emergence)  its 
intex’action  with  its  surroundings  is  held  to  be  absolutely  vital.  In  fact 
the  buoyancy  force  is  balanced  by  form  drag  (produced  by  flow  around  the 
bubble  cap)  and  is  not  producing  a resultant  upward  acceleration.  On  the 
contrary,  the  bubble's  upward  velocity  is  in  general  decreasing  due  to  ero- 
sion. To  summarize:  the  parcel  method  prescribes  no  relations  between 
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the  velocity  and  radius  of  a convective  element  and  an  erroneous  relation 
between  velocity,  buoyancy  and  time.  The  present  model  suggests  quantita- 
tive relations  between  these  parameters,  which  are  tested  here. 

Although  the  wake  production  process  has  not  been  explicitly  examined, 
it  is  hypothesized  to  be  the  mechanism  by  which  mixing  with  the  surroundings 
and  entrainment  of  outside  air  into  the  main  cloud  body  are  occurring.  This 
topic  must  be  the  subject  of  much  future  work,  once  the  simpler  features  of 
the  present  model  nave  been  established  quantitatively.  Meanwhile  certain 
further  deductions  concerning  interaction  of  bubbles  and  environment  may  be 
drawn  from  the  uresent  data. 

c)  The  Effects  of  Wind  Shear  and  the  Horizontal  Drag  on  a Bubble 

In  the  case  of  a bubble  ascending  through  a shearing  wind  field  it  was 
shown  in  I,  Appendix  II,  that  actually  equation  (5)  should  be  replaced  by 

3 JL  _ (-t)^  gB  (1  - w/gB)t 

~ E ~ = G = (6) 

2 v + ,,2^ 

where  u is  the  horizontal  velocity  of  the  bubble  relative  to  the  air.  Tt  was 
further  shown  that  in  most  shearing  fields  the  difference  between  the  re- 
sults of  (£)  and  of  (6)  would  be  only  about  2%,  but  that  when  u is  of  the 
same  size  as  w,  use  of  (5)  gives  about  a 20*  overestimation  of  G.  Tables  I 
and  III  indicate  that  in  the  present  case  u is  indeed  comparable  to  w,  since 
it  ranges  from  9-1$  raps  in  the  region  of  strong  shear.  However,  it  was  also 
shown  that  a rate  of  buoyancy  decrease  of  the  magnitude  experienced  by  the 
present  bubbles  at  these  heights  should  lead  to  about  a 1^20%  underestima- 
tion of  G.  Since  these  errors  are  compensatory  and  since  the  present 
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observations  are  probably  not  reliable  to  better  than  15-20%  ary  how,  no  at- 
bdUp  \j  will  be  made  here  to  introduce  such  refinements. 

Sines  vc  have  in  Tables  I and  III  and  in  Figure  6 some  measurements  of 
the  horizontal  tower  velocity,  u,  and  also  the  external  wind  profile,  io  may 
be  possible  also  to  test  a horizontal  equation  of  motion  for  the  cloud  tow- 
ers. In  I,  Appendix  H,  is  suggested  an  equation  analogous  to  the  vertical 
equation  of  motion  with  wind  shear  present,  namely 


% ~ ~ ,f~2 21 

/u  + v 


K (u2  + w2) 


where  Uu  is  the  component  of  the  bubble's  horizontal  motion  alouK  the  axis 

U 

determined  by  the  external  wind,  Ugj  u is  tne  relative  motion  between  bubble 
ana  surroundings,  namely  u = _ Ugj  and  the  dot  indicates  time  differen- 

tiation. The  constant  K was  shown  to  be  the  same  as  K in  equation  (3), 
namely  9/Ufi.  Subtracting  Ug  from  both  sides  of  (7),  we  have 


u = - wU'  - ~r  2 " 2>  K (u2  + w2) 

vn  + w£ 


since  Ug  = wU'  = wdu/dz  and  u = - Ug  and  U'  or  du/dz  is  the  external  wind 

shear. 

If  the  bubble  has  been  rising  through  constant  shear  for  several  hundred 
meters  and  if  it  still  is  fairly  near  the  beginning  of  its  life  cycle  so  that 
K is  changing  only  slowly,  u is  very  small,  as  is  w,  and  (8)  may  be  solved 
approximately  to  give 


9 y u 2 -r  w* 


after  substitution  of  (3a)  for  K. 

The  region  of  strong  shear  on  Figure  6 corresponds  to  the  earlier  phases 

of  the  cycles  o f the  tower-  -studied.  Hence  (9)  may  be  tested.  Table  I shows 

that  Tower  X possessed  a u of  about  -lf>  raps  in  this  region  and  Table  III 

gives  an  average  for  Tower  Q of  about  -12  mps.  Assuming  that  these  bubbles 

had  been  rising  through  the  shearing  layer  long  enough  lor  u to  be  very  small, 

-2  _1 

and  taking  the  observed  U'  of  3 x 10  sec  , R for  Tower  X comes  out  1.8  km 
and  R for  Tower  Q comes  out  l.U  km,  compared  to  the  calculated  values  of  2.1 
and  1.3  km  respectively  which  appeared  in  Table  VII.  These  independent  re- 
sults are  consistent  with  the  predictions  of  the  model  and  agree  within  ob- 
servational error,  although  it  must  bs  conceded  that  the  measurements  and 
the  approximations  made  are  such  that  this  high  degree  of  agreement  is  some- 
what surprising. 

It  is  not  clear  from  the  observations  whether,  as  the  bubble  continues 
its  life  cycle,  the  horizontal  drag  coefficient  increases  in  the  same  manner 
as  its  vertical  counterpart,  which  the  model  declares  that  it  should,  although 
this  is  at  least  roughly  suggested  by  the  fact  that  u = - ug  seems  to  be 

decreasing  in  magnitude  before  evaporation  sets  in  and  makes  it  impossible  to 
measure  further.  Changing  wind  shear  coupled  with  the  difficulty  of  tracing 
a single  point  in  its  horizontal  travel  precludes  an  accurate  quantitative 
test  of  (8)  solved  for  E (after  substitution  of  K = 9/lfR  and  R = - EgBt) 
such  as  was  done  for  the  vertical  equation  (3). 

In  a shearing  wind  field  a rising  bubble  should  have  a horizontal  as 
well  as  a vertical  component  to  its  trailing  turbulent  wake.  Earlier  obser- 
vations by  Malkus  (1949)  suggest  the  presence  of  such  a wake  region  in  the 
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horizontal,  since  several  hundred  feet  of  "dying  cloud*  was  commonly  found  on 
the  downshear  edge  of  cumuli. 


V.  CONCLUSION 

The  clouds  studied  in  this  paper  are  of  interest  for  two  main  reasons: 
first,  they  are  the  largest  clouds  in  conditions  of  strongest  shear  to  which 
the  techniques  described  here  have  been  applied.  Previous  clouds  investi- 
gated by  this  project  have  been  small  or  middle-sized  trade  cumulus  and  fair 
weather  cumulus  of  middle  latitudes,  none  of  whose  tops  reached  the  freezing 
level.  Some  similarities  and  some  contrasts  between  these  clouds  and  the 
smaller  ones  studied  earlier  are  striking.  The  outstanding  similarity  is 
the  applicability  of  the  same  erosion  constant,  E,  tc-  the  towers.  The  pre- 
dominant contrast  is  the  relatively  less  important  role  played  by  mixing 
with  the  clear-air  environment  in  the  bigger  clouds.  This  is  not  to  say 
that  interaction  between  cloud  elements  and  their  surroundings  is  any  less 
important,  since  it  was  demonstrated  here  that  the  strongest  bubbles  studied 
emerged  from  cloud  top  with  velocities  less  than  1/6  that  of  a non-interact- 
ing parcel.  In  fact,  it,  has  been  indicated  that  a minimum  cloud  dimension 
is  necessary  to  shield  the  innermost  core  from  dilution.  In  the  case  of 
these  clouds,  a mass  roughly  9 km  across  (in  the  plane  of  the  wind  and  the 
wind  shear)  could  produce  at  intervals  of  about  20  minutes,  undiluted  bub- 
bles 2 km  in  diameter  which  emerged  from  cloud  top  rising  at  11  mps.  The 
relation  between  radius,  buoyancy,  velocity,  and  lifetime  of  these  elements 
agreed  very  satisfactorily  with  the  relations  established  by  the  work  of 
Malkus  and  Scorer  in  paper  1. 

The  greatest  importance  meteorologically  of  the  clouds  studied  is. 
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however,  their  role  in  the  over-all  picture.  It  must  be  supposed  that  clouds 
of  this  type  and  scale,  although  rarely  observed  from  the  ground  due  to  pre- 
vailing overcasts  under  disturbed  conditions,  are  the  rule  and  not  the  excep- 
tion when  an  extended  polar  trough  reaches  down  to  trade-wind  latitudes. 

It  is  well  established  (Riehl,  1950)  that  the  strong  and  high-penetrating 
convective  activity  which  occurs  in  these  troughs  is  the  means  by  which  ths 
upper  troposphere  receives  the  major  fraction  of  its  latent  energy  supply. 

A sample  calculation  will  illustrate  the  magnitudes  involved.  Assume  that  at 
the  25,000  ft.,  level,  a single  saturated  cloud  tower  is  active  at  a given 
moment  in  an  area  30  km  on  a side  (roughly  correct  for  April  1,  1953).  If 
the  updraft  speed  is  IQ  zips  and  the  dr* it  is  2 km  on  a side,  continuity  will 


be  met  by  subsidence  at  li  cm/sec  over  the  remaining  area.  Since  even  the 
more  dilute  towers  at  this  elevation  have  about  3.0  gm/kg  water  vapor  con- 
tent and  the  ambient  air  from  the  sounding  less  than  0.5  gm/kg,  the  net  up- 
ward transport  of  water  vapor  is  probably  not  less  than  6 x 10^  gm/sec  or 
3.6  x 10^  cal/sec  in  latent  heat.  Distributed  over  the  area  involved,  this 
is  U00  cal/cm  per  day.  This  figure  should  be  compared  to  a several  months’ 
average  evaporation  rate  of  250  cal/cnf  per  day  calculated  lor  a similar 
latitude  in  the  Pacific  trade-wind  region  by  Riehl,  Teh,  Malkus,  and  La  Seur 
(1951).  Clouds  such  as  these  thus  play  a major  role  in  maintenance  of  the 
energy  budget  of  the  trades  and  on  that  basis  alone  merit  further  stuejy. 
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up  this  particular  section  of  the  data,  and  the  criticism  and  valuable  ideas 
provided  by  Dr.  R.  S.  Scorer  were  essential  in  the  analysis. 
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Append^ x 

Concerning  the  Applicability  of  the  San  Juan  Radiosonde  and 
Radio -Wind  Data  to  the  Location  of  the  Clouds  Studied 

Since  the  clouds  are  160  miles  away  from  San  Juan,  and  the  midpoint 
of  the  observations  upon  them  occurred  two  hours  subsequent  to  the  1500 
GOT  sounding,  the  important  question  arises  whether  the  wind  and  tempera- 
ture structure  at  the  time  and  place  of  the  clouds  might  not  differ  sig- 
nificantly from  the  San  Juan  sounding  which  was  used  in  the  calculations. 

Since  it  can  be  established  by  investigation  of  the  synoptic  situa- 
tion on  this  and  previous  days,  that  local  changes  occurring  in  these  two 
hours  are  probably  insignificant,  the  question  concerns  the  magnitude  of 
advective  changes. 

By  using  the  mean  wind  velocity  throughout  the  levels  of  interest 
(9.U-lii  km)  and  projecting  backward,  the  air  in  the  cloud  locality  at  the 
middle  of  the  run  was  100  miles  due  north  of  San  Juan  at  the  time  of  the 
radiosonde.  If,  for  example,  the  weak  stable  layer  occurring  between 
about  10-11  km  (which  coincides  with  the  layer  of  high  wind  shear)  had  a 
slope  northward  of  l/l!>0,  or  that  common  for  a middle-latitude  frontal 
system,  it  would  be  at  a 3500  ft.  or  approximately  10$  different  elevation 
100  miles  to  the  north,  or*  in  other  words,  the  sounding  at  the  clouds 
would  be  shifted  in  height  by  10$  at  the  time  and  place  of  the  observa- 
tions upon  ther . In  the  trade-wind  region,  a north-south  slope  this  great 
■is  almost  surely  out  of  the  question. 

The  results  of  all  the  calculations  taken  together  further  support 
the  contention  in  this  Appendix.  A .'.0$  error  in  sounding  heights  relative 
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to  cloud  heights  would  give  buoyancy  errors  of  about  20%  in  the  region  from 
about  9.7-U.O  km  where  the  buoyancy  is  decreasing  rapidly  with  height,  and 
no  error  from  there  up  to  near  the  level  where  the  bubbles  ended  their  cy- 
cles. The  constancy  of  u from  9.7  km  on  up  precludes  the  possibility  of 
such  an  error,  either  in  sounding,  or  in  cloud  heights,  unless  by  chance 
they  should  both  be  in  error  in  a compensating  manner.  The  additional  fact 
that  the  wind  sounding  began  to  shear  at  just  the  height  that  the  clouds 
began  to  show  shear  on  the  film  further  mitigates  against  any  but  a fortui- 
tously compensatory  error. 
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Titles  for  Illustrations 

Fig.  1.  Time  section  for  San  Juan,  Puerto  Rico,  from  1200  GCT  March  31, 

1953  to  0600  GOT  April  2,  1953.  The  ordinate  is  height  in  feet.  Winds 
are  plotted  in  knots,  a short  barb  indicating  5 knots,  a long  barb  10 
knots,  and  a solid  triangle  50  knots.  The  solid  lines  arc  isotachs  of 
wind  speed,  drawn  every  10  knots,  with  J indicating  maxima.  The  dashed 
line  indicates  the  base  of  the  westerlies,  and  the  dotted  line  shows  the 
top  of  the  moist  layer.  The  heavy  solid  line  marked  PT  indicates  the 
polar  trough  which  passed  San  Juan  moving  eastward  on  March  31,  and  the 
solid  lane  marked  EW  denotes  a weak  easterly  wave  which  passed  San  Juan 
moving  westward  early  on  April  1. 

Fig.  2.  Still  photograph  from  the  PBY  showing  the  clouds  studied.  The  ar- 
rows indicate  the  points  upon  which  triangulation  was  performed.  The 
aircraft  was  being  flown  at  60U0  ft.  toward  293°  True  and  its  location 
is  given  in  Figure  5. 

Fig.  3.  Selected  prints  from  the  motion  picture  camera  on  Anegada  from 
which  the  calculations  were  made.  The  bubble  followed  on  Tower  Y is 
shown  by  the  line;  Tower  W is  indicated  by  a square;  Tower  C by  the 
arrow  alone;  Tower  Q by  the  inverted  triangle;  the  streamer  on  Cloud 
II  by  the  erect  triangle,  and  Tower  X by  the  X.  The  camera  was  pointed 
toward  22°  T throughout  the  run.  The  print  number  is  denoted  beneath 
each  frame  and  consecutive  prints  were  1.2  minutes  apart  (every  twenty 
xi  cuuco  ui  uhc  mot  ion  picture  film,  which  exposed  at  the  rate  of 


one  frame  every  3.6  seconds) 
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FLg.  U.  Additional  selected  prints  from  the  motion  picture  camera  on  Anegada* 
On  Print  1*0,  Tower  P is  shown  by  the  arrow  alone;  Tower  X is  indicated 
by  an  X,  and  the  streamer  on  Cloud  H by  an  erect  triangle.  Print  hli 
was  the  one  used  in  triangulation  (compare  to  FLg.  2).  The  arrows  alone 
indicate  the  points  triangulated  upon  (see  results  in  Pig.  $).  The  ar- 
row alone  on  Prints  $2  and  58  (the  last  one  in  the  run)  shows  Tower  R. 

The  dissipation  of  the  entire  cloud  mass  may  be  noted  on  the  last  print. 

Fig.  5*  Scaled  map  showing  results  of  the  triangulation  using  Figure  2 and 
Print  1*1*,  Figure  1*.  The  camera  on  Anegada  was  located  at  A,  pointing 
in  the  direction  of  the  dashed  arrow  (toward  22°).  The  FBI  was  located 
at  B,  heading  in  the  direction  indicated  by  the  solid  arrow.  The  three 
points  indicated  by  arrows  on  Cloud  T (Fig.  2 and  Fig.  1*)  lie  along  the 
line  marked  "Cloud  I"  which  lies  very  closely  in  the  plane  290o-110° 
and  78  miles  perpendicular  distance  from  A.  The  single  point  triangu- 
lated upon  in  Cloud  I lies  66  miles  perpendicularly  from  the  plane 
290°-110°  through  A. 

Fig.  6.  The  wind  profile  plotted  from  the  San  Juan  radio-wind  observation, 
1600  GCT,  April  1,  1953.  These  points,  marked  by  circles,  were  ob- 
tained upon  reduction  of  the  original  data  for  this  purpose.  The  curve 
on  the  left  is  windspeed  in  meters  per  second,  that  on  the  right  is 
direction  in  degrees  from  true  North.  The  layer  of  strong  shear  has 
been  carefully  checked  and  found  on  the  other  soundings  in  the  area, 
and  on  previous  and  subsequent  soundings  at  San  Juan.  The  other 
points,  triangles,  bars,  squares,  and  X's  are  horizontal  speeds  of  the 
various  cloud  towers,  as  denoted  on  the  figure. 
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Fig.  7.  Tephigram  showing  San  Juan  radiosonde  observation,  1500  GOT,  April 
1 , 19$ 3,  reduced  from  the  original  data  to  obtain  more  closely  spaced 
points  than  that  of  the  teletype  transmission.  The  sounding  is  given 
by  the  solid  line.  The  dashed  line  is  the  moist  adiabat  intersecting 
the  sounding  at  cloud  base  (9^0  mb)  used  to  obtain  the  buoyancy  for 
the  strongest  bubbles,  and  the  dotted  line  is  the  moist  adiabat  used  in 
obtaining  the  buoyancy  of  some  of  the  weaker  bubbles. 

Fig.  8.  Height  of  top  of  bubble  in  Tower  X in  km  as  a function  of  time  in 
seconds.  The  time,  -t,  is  chosen  as  zero  when  the  bubble  levels  off 
and  becomes  fuzzy  in  appearance  on  the  film.  Several  points  commonly 
appear  on  the  same  horizontal  because  readings  could  only  be  made  to 
the  nearest  0.0 $ cm  on  the  microfilm  viewer.  Readings  were  made  both 
in  a forward  and  backward  direction  on  the  film  and  were  found  to 
coincide  to  this  degree  of  accuracy.  Velocities  were  obtained  by  read- 
ing first  differences  of  the  heights,  and  accelerations  by  taking  sec- 
ond differences  (see  inset  graph).  Accelerations  were  checked  by  the 
changes  in  velocity.  The  calculations  made  from  these  curves  are  pre- 
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